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Abstract 

>, 

Photons may convert into axion like particles in the magnetic field of active galaxies, 

00 . 

. traverse the Universe and then convert back into photons in the magnetic field of the 

Milky Way. This is a potential explanation for the candidate neutral ultra-high-energy 
{E > 10 18 eV) particles from distant BL Lac type objects which have been observed by 
the High Resolution Fly's Eye experiment. We find that there is a correlation between 
the direction of these sources and the regions on the sky where the photon-axion mixing 



probability is highest. The probability of obtaining this correlation by chance is only 
2.4%. 



1 Introduction 

Axions are pseudo-scalar particles which arise as the Nambu-Goldstone Bosons of the broken 
Peccei-Quinn symmetry [lj. They obtain a mass when the CP violating QCD theta term 
is driven to zero in agreement with observations [21 |3j. When motivated in this way, the 
relationship between the axion mass and coupling is related to the pion mass and decay 
constant such that for a given axion mass, the coupling to photons is determined up to 
factors of order a few [mM ~ vn^f^ where m is the axion mass and M is the inverse 
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axion coupling, see section [2]). While considerable experimental and theoretical work has 
eliminated much of the parameter space of such models, axions are still a viable candidate 
for both the solution of the strong-CP problem and for cold dark matter. 

The term Axion Like Particle refers to a particle with a similar Lagrangian structure 
to the Peccei-Quinn axion but where the constraints on the parameters of the Lagrangian 
have been relaxed. In other words there may be particles like the axion weakly coupled to 
the standard model even if they do not solve the strong-CP problem. For simplicity, in the 
rest of this paper, we shall refer to all such particles as axions, while the particular kind 
of particle associated with the solution of the strong-CP problem we shall refer to as the 
Peccei-Quinn axion. 

Axions have been invoked to solve a variety of different problems in physics and astro- 
physics. For example, it has suggested that they might be responsible for the dimming of 
supernovae - photons from distant type la supernovae might convert into axions as they 
cross the Universe to reach us which may explain the apparent low luminosity of high red- 
shift supernovae normally subscribed to the presence of a cosmological constant [3]. Such 
models are interesting, although there may be problems with the frequency dependence of 
the dimming effect and they appear difficult to reconcile with baryon acoustic oscillation 
observations [3, El [7] . 

Axions were also invoked to explain the anomalous results from the PVLAS collaboration 
which suggested possible optical activity of the vacuum in the presence of magnetic fields 
[S]. Subsequent results from the PVLAS experiment retracted these initial claims however, 
negating the need for such a particle [9]. 

In Ref. [10J, photon-axion oscillations in intergalactic space have been suggested as an 
explanation of super-GZK cosmic rays detected by the AGASA experiment; mixing in the 
source was not considered. Since such mixing means that photons spend some of their time 
as axions while on route to earth, the attenuation length of photons is effectively increased. 
Unfortunately it seems that for the parameters of Ref. [10], the original flux of photons in 
the source should exceed the flux observed at the Earth by several orders of magnitude - all 
these additional photons have to lose their energy in cascades on the background radiation 
which results in a conflict with EGRET limits on the diffuse gamma-ray background. 

More recently, the detection of TeV photons from objects at cosmological distances has 
led to a new reconsideration of axions. It is difficult to explain how such photons could 
reach earth given the opacity of the Universe at those wavelengths due to pair production on 
the background infrared radiation. It has been suggested that the mixing of photons with 
axions in the intergalactic magnetic field may explain this, although the required intergalactic 
magnetic field has to be on the high side [IT] . Another suggestion is that photons are 
converted into axions in the magnetic field of the active galaxy itself, which is a rather 
reasonable assumption for axions with low masses. If such a mixing were to take place 
efficiently, up to one third of the initial high energy photon flux may cross the Universe in 
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the form of axions before being converted back into photons in the magnetic field of the 
Milky Way, avoiding the attenuation that photons would experience as they travel across 
the Universe. Since the authors of [12] identified the axion parameters and galactic magnetic 
field which can explain the arrival of TeV photons from cosmological sources, we shall test 
the same parameters and magnetic field to see if they might also explain the origin of these 
apparently neutrally charged ultra high energy cosmic rays. 

One of the most fascinating predictions of theories which contain axions is the idea that 
one may 'shine light through walls' by converting photons to and from axions on either side 
of a wall using strategically placed magnetic fields. In this work we are doing the same 
experiment but we are using the Universe as our wall and galaxies for our magnetic fields. 

In section [2] we shall go over the mathematics of the mixing phenomenon and discuss the 
mixing of photons and axions in astrophysical sources. Then in section [3] we will discuss the 
evidence for the arrival of ultra high energy cosmic rays from directions coincident with BL 
Lac objects before discussing photon-axion mixing as a possible explanation for these events 
in section HI We will then discuss the magnetic field of the Milky Way that we will assume 
to perform our analysis in section [5] before preforming the search for correlations between 
regions of large B-field and cosmic ray events in section [6j Finally we will list some of the 
consequences of this model and other ways to test it before moving to our conclusions. 

2 Photon-axion mixing in astrophysical objects. 

The Lagrangian describing the photon and axion takes the following form (similar results 
hold for a scalar), 

£ = -(d^ad.a - mV) - \^ 7 F lll/ F^ - -F^F^, 
2 v m 4 M 4 

where F pv is the electromagnetic stress tensor and F pv = e pup \F p \ is its dual, a denotes 
the pseudo-scalar axion, m is the axion mass and M is the inverse axion-photon coupling. 
Because of the F I1V F^ U term, there is a finite probability for the photon to mix with the axion 
in the presence of a magnetic field. Mixing also occurs between photon components with 
different polarizations [T3T fT5"| ITS]. We will be interested in light, m < 10~ 7 eV, axions with 
inverse coupling mass scale M ~ fewxl0 10 GeV which are the parameters that can explain 
the observations of TeV photons from cosmological sources [12] . For axions of these masses 
the most stringent bound on the coupling, M > 1.1 x 10 10 GeV at the 95% CL, has been 
placed by the CAST experiment [T7] . 

Technically, the mixing may be described as follows. We represent the photon field 
A(t,x) as a superposition of fixed-energy components A(x)e~ tU}t . If the magnetic field does 
not change significantly on the photon wavelength scale and the index of refraction of the 
medium |n — 1| <C 1, one can decompose [16] the operators in the field equations as (for a 
photon moving in the z direction) uo 2 + d 2 — *> 2uo{uo — id z ), so that the field equations become 
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the mixing is determined by the refraction parameter A p , axion-mass parameter A m and 
mixing parameter Am, which are equal to 
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respectively. Here <J^ = 4iran e /m e is the plasma frequency squared (effective photon mass 
squared), n e is the electron density, B i: i — x,y are the the components of the magnetic 
field B at right angles to the direction of propagation, m e is the electron mass, a is the 
fine-structure constant and u is the photon (axion) energy. For constant magnetic field and 
electron density, the conversion probability is 
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where 



A L = (A P - A m f + 4A M . 



If B and n e change spatially, the probability can be found by a numerical solution of Eqns. ([T]). 

Looking at Eq. one sees that the conversion probability is large in the maximal- mixing 
regime when 

(A p - A m f < 4A M . 

In real astrophysical sources, the precise cancellation of A p and A m will occur only in small 
regions, so the maximal- mixing conditions, which both must be met, are 
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which are equivalent to 
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Figure 1 : Typical values of the magnetic field [H] and electron density [T9J, EH] in various 
astrophysical objects (IG: intergalactic space, MW: the Milky Way, NLR and BLR: narrow- 
and broad-line regions in active galactic nuclei). The condition (Jlj) is satisfied above the 
thick line for energies u > 10 19 eV, above the dashed line for u > 1 TeV and above the thin 
line for uj > 10 MeV (for M = 10 10 GeV). 



-<< 1020 («)(§)• < 4 » 

In addition, to have large mixing one should require that the size L of the region in which 
conditions (EJ), (HI) are fulfilled should exceed the oscillation length, 



that is 

£ > 5 . 8X10 -3 PC (§H_J1_). (5) 

From Fig. [TJ one sees that Eq. (Jlj) is certainly fulfilled for ultra-high-energy particles in 
all astrophysical gamma-ray sources. For axion-photon coupling close to its experimental 
limit, the condition (jlj) is met down to energies as low as ~ 10 MeV. The condition ([3]) is 
illustrated in Fig. El The condition (jSj) is the most restrictive. It depends on both the size 
and the magnitude of the magnetic field and can be superimposed [21] on the Hillas plot 
for various astrophysical sources (Fig. [3]). We see that if an axion like particle exists with 
the mass and coupling outlined above, high-energy photons intensively mix with it in many 
gamma-ray sources. As a result, the axion flux F a = F^/2 accompanies the gamma-ray flux 
F 1 independently of the gamma-ray emission mechanism (for the maximal mixing, fluxes of 
axions and of photons of each polarisation are equal). 
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Figure 2: The condition ([3]) on the parameter plane "photon energy" - "magnetic field". 
The condition is satisfied above red lines (dotted: for m = 10~ 7 eV, thick for m = 10~ 9 eV, 
thin for m = 10 -11 eV). Horizontal lines indicate typical B values for various astrophysical 
sources, cf. Fig. [T]). Thick blue (yellow for the Milky Way) parts of the horizontal lines 
correspond to the energy range where the condition is satisfied for m = 10~ 9 eV (for 
M = 10 10 GeV). 
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Figure 3: The condition (jSJ) for M = 10 10 on the updated Hillas plot [18J. The condition 
is satisfied in the shadowed region. The Milky Way parameters are denoted by a star. Also 
shown are parameters for anomalous X-ray pulsars and magnetars (AXP), neutron stars 
(NS), central black holes (BH) and central parsecs (AD) of active galaxies (low-power Seyfert 
galaxies (Sy) and powerful radio galaxies (RG) and blazars (BL)), relativistic jets, knots (K), 
hot spots (HS) and lobes (L) of powerful active galaxies (RG and BL); non-relativistic jets of 
low-power galaxies (Sy); starburst galaxies; gamma-ray bursts (GRB); galaxy clusters and 
intercluster voids. 
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3 Neutral particles from distant sources 



A number of studies suggest that a correlation may exist between the arrival directions 
of cosmic rays and catalogues of BL Lac objects. This correlation exists without taking 
into account the magnetic field of the galaxy, suggesting that the cosmic rays experience 
zero deflection as they traverse this field and are therefore neutral particles, challenging 
conventional models of cosmic-ray physics. 

These claims are based upon two samples of cosmic rays, the first sample combines events 
from the Akeno Giant Air-Shower Array (AGASA) of cosmic rays with estimated primary 
energies E > 4.8 x 10 19 eV) and a sample from the Yakutsk Extensive Air Shower Array 
(Yakutsk) of events with estimate primary energy E > 2.4 x 10 19 eV. An excess of correlations 
between the position of BL Lacs and the arrival direction of cosmic rays in this combined 
data set was seen at separations less than 2.5° |22j . 

Similarly, a sample of events with E > 10 19 eV observed by the High Resolution Fly's 
Eye detector (HiRes) tested positive for correlations between source and BL Lac objects at 
angular separations less than 0.8° [23J. 

In both cases the separation was consistent with the detector's angular resolution (which 
was much better in HiRes than in AGASA and Yakutsk). The correlation with the HiRes 
sample was confirmed in an unbinned study and was found to extend to lower energies |24J. 
The probability to observe the correlation with three independent experiments by chance 
was estimated by [25] as 3 x 10~ 5 by a Monte-Carlo study. 

The correlation between BL Lacs and UHECRs seen in HiRes data [23J has been tested 
recently by the Pierre Auger Collaboration [2E] arid no positive signal has yet been found. 
However, it turns out that this is unsurprising for the following three reasons:- 

Firstly, Pierre Auger is located in the Southern hemisphere and sees different BL Lacs to 
other experiments and due to incompleteness of the astronomical catalogs, fewer potential 
UHECR emitters are known in the South. 

Secondly the angular resolution of the Pierre Auger array is also inferior to that of the 
Hi- Res array which means that the sensitivity to such correlations can only be achieved with 
much more data. 

Finally, as has been pointed out in [23] and [24] and further discussed in [27J, the correla- 
tion observed by HiRes implies neutral cosmic particles traveling for cosmological distances, 
a fact which requires unconventional physics. Most probably the primary particles of the 
resulting air showers are neither protons, nor nuclei. However, the energy determination of 
the PA surface detector is extremely sensitive to the type of the primary cosmic particle 
because of very strong sensitivity of water tanks to muons in the air shower. For instance, 
energies of gamma rays are always underestimated by a factor of a few (see e.g. [2E1I2H])- Due 
to the steeply falling spectrum of UHECRs, this may dilute the observed signal. It would 
therefore be interesting if the Pierre Auger collaboration were able search for the correlation 
using their Fluorescence detectors rather than the water tanks. 
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Having discussed the evidence for the correlation between the arrival direction of ultra 
high energy cosmic rays and BL Lac objects, we will move on to look at the use of axions to 
explain how neutral particles could traverse the Universe without complete attenuation. 

4 Axions as ultra high energy cosmic rays. 

In the framework of the Standard model of particle physics and assuming standard astro- 
physics, neutral particles with energies ~ 10 18 eV cannot propagate for > 100 Mpc, the 
distance to the nearest BL Lacs. The only exception is neutrino which can be excluded as 
an explanation for these events by considering the height of development of the atmospheric 
showers and noting that they are not close enough to the ground to be consistent with the 
weak interaction cross sections. 

Photons interact with the background radiation which results in pair production and the 
development of electromagnetic cascades. Known unstable particles decay at much shorter 
distances (Fig. HI). In the framework of more involved descriptions which do not require new 
physics, the neutral particles may be created in interactions of protons inside or not far from 
the Milky Way; however, in this case the observed effect also cannot be explained [27J. 

Even beyond the Standard model it is difficult to find a non-contradictory explanation of 
the observed correlations. New stable strongly interacting particles [2B E2] should be heavy 
enough not to be detected in accelerators, but the probability to create such a heavy particle 
is low and therefore one would expect for each one of these particles that there should be 
huge fluxes of accompanying radiation in conflict with constraints on diffuse cosmic gamma 
radiation. 

Models which suggest the existence of a relatively heavy (~MeV) axion like particle 
(sgoldstino) [33] suffer the same problem. 

In the models where there is an enhanced neutrino-air cross section (see e.g. [21]), besides 
some theoretical difficulties, the cross section rise is not sufficient to explain the shower 
development. Only in the models with Lorentz-invariance violation [35] decaying neutral 
particles (neutron or 7r° meson) might be stable in a certain energy region and propagate to 
cosmological distances. It might be argued however that postulating the existence of a new 
particle is less drastic than altering the framework of relativity. 

The scenario which we investigate here is based on the mixing of photon with light 
axions. The parameters of the model which work in explaining the conundrum of the neutral 
primaries outlined above do not contradict with any experimental limits and may allow one 
to explain some other astrophysical puzzles as well. 

The maximal mixing conditions ©, ©, © satisfy (cf. Figs. [Q [21 E]) for m < 10~ 7 eV 
and M > 10 10 GeV both for BL Lacs and for the Milky Way. At the same time, for the 
magnetic-field values ~ 10 -11 G in inter-cluster voids [36], the condition ([5]) is not satisfied, 
so the following scenario is possible:- In the source, protons are accelerated to ultra-high 
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log 10 (E/eV) 

Figure 4: Attenuation length of different kinds of particles as compared to the distance 
to the nearest BL Lac (blue horizontal line) and the size of the Universe (upper bound 
of the plot). The green line corresponds to photons, blue to neutrons and red (shown for 
comparison) to protons. The lines for protons and photons are taken from the review |30j. 
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(E > 10 20 eV) energies (according to Ref. [18] . acceleration of protons in BL Lacs up to 
these energies contradicts neither the Hillas criterion nor the radiation losses). Interaction 
of these protons with the intense blazar emission results in the pion photo-production similar 
to the GZK effect which for a fraction of the accelerated particles takes place directly in the 
source. Intensive mixing in the source magnetic field converts 1/3 of photons into the axions 
of the same energy so the source emits the axion-photon beam into space towards earth. 
Further mixing in intergalactic space before the photon-axion beam arrives at the Milky 
Way is suppressed due to the small magnetic fields in that region. The photon part of the 
beam interacts with background photons and loses energy while the axion part propagates 
unattenuated. Then, upon arrival at the Milky Way where the magnetic field is several 
orders of magnitude higher than the intergalactic field, intensive mixing again takes place 
and a significant fraction (2/3 for maximal mixing) of the axions are converted back into 
photons which are then detected as neutral particles from BL Lacs. The maximum fraction 
of photons detected in cosmic ray detectors on earth can be 2/9 of the total flux of photons 
of same energy emitted in the source. 

Because the magnetic field of the Milky Way is non-uniform, we expect the probability 
for axions which arrive from far away to convert back into photons to depend upon the 
direction they move through the galactic magnetic field before the arrive at Earth. This is 
the subject of the next section. 

5 Direction dependence of mixing in the Milky Way 

Magnetic fields in the Milky Way typically have [iG magnitudes and the structure of the 
fields are thought to be strongly related to the structure of the spiral arms in the galaxy. It is 
difficult to reconstruct the detailed field strength and structure from the present observations, 
which are mostly based on the Faraday rotation measurements. The exact structure of the 
field is still uncertain, although there are several phenomenological models proposed in the 
literature [33 M, [39] . 

The same effect that we are presenting in this paper as a possible explanation for the 
origin of neutral ultra high energy cosmic rays can also explain the observation of TeV 
photons from Blazars, photons which should not be able to cross cosmological distances due 
to their interaction with the cosmic infrared radiation [12]. Therefore, for consistency with 
|12j . we concentrate on the magnetic field model which allows one to explain the arrival of 
TeV photons. This turns out to be the Harari, Mollerach and Roulet (HMR) model [37]. It 
is expected that the other two models would give worse agreement with the arrival directions 
of the observed events, as it was found in Ref [12J. 

The HMR model is described by the following equations 

B e = B(r,6,z)cos(p) , B r = B(r, 9, z) sin(p) . (6) 
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Figure 5: Regions of high probability of mixing for high energy photons passing through 
the HMR field described in the text as a function of position in the sky. 



where (6, r, z) are cylindrical Galactic (Galactocentric) coordinates. The strength of the field 
is given by 

B(r,e,z) = — tanh 3 ( — J cos [9 — In f ) ] x 

r \PiJ V tanp \£ J J 

X (2cosh(z/ 2l ) + 2cosh(^/z 2 )) ^ G (7) 
The free parameters are taken to be 

Z\ = 0.3 kpc , z 2 = 4 kpc , p = —10° , £ = 10.55 kpc , p\ = 2 kpc . (8) 

Such that p represents the angle between the opening angle of the spiral arms. The distance 
from the Sun to the Galactic center is R = 8.5 kpc. We also assume that the magnetic field 
is zero at radii r > 20kpc. 

We then evolve the photon-axion system using equation ([T]). In order to achieve this we 
have to move to galactic solar-centric coordinates (r 3 , 1, b) using the appropriate transforma- 
tions where I and b are the direction on the sky and r 3 is the distance along that direction. 
We then obtain the integrated probability as a function of (I, b) which is plotted in figure [5] 

In [12] the authors do not evolve the full equations but rather argue that a good enough 
approximation given that the maximal mixing criteria are fulfilled is to evaluate | j B±dr 3 \ 2 
along the appropriate direction. We perform both the full method and the method used 
in that paper and are happy to confirm that their approximation seems to be justified, 
nevertheless we use the more complete method in our work, since we have the code to do it. 
We also seem to agree with the magnitude of the mixing probability in the aforementioned 
work at least to within factors of order one but since there is no Lagrangian in that paper 
we are unsure precisely how their parameter gn relates to our 1/M. 

A final issue worth mentioning is the fact that we do not include turbulence in the mag- 
netic field. We predict that this should not be a problem since the amplitude of the turbulent 
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events vs. exposure events vs. exposure x mixing 

Figure 6: In the figure on the left we plot the exposure of the HiRes cosmic ray detector - 
darker shaded regions have had more exposure. The red data points correspond to cosmic 
rays correlated with the position of BL Lacs. On the right we plot the same events on the 
background of the combined probability of exposure multiplies by mixing probability due to 
an axion and magnetic field of the form specified in the text. 

field should follow on average the background uniform field. There is the possibility that 
short scale fluctuations could suppress the mixing slightly because the mixing length is com- 
parable to the sizes over which such turbulence will occur. However the lack of suppression 
in mixing between our fuller treatment relative to the previous averaged treatment of [T2] 
encourages us to believe that this will not be a huge effect. One result of having a turbulent 
field may be to increase the photon-axion mixing probability in the Galactic plane which 
might change the correlation. It may be worth checking this in more detail in the future. 

Using the results of this section and those of section [2] we have therefore been able to 
obtain a probability of axion-photon mixing which is a function of which direction we look 
through the magnetic field of the Milky Way. We now need to find out if the arrival directions 
of cosmic rays are correlated with those directions in which there is a large probability for 
mixing. 

6 Results 

For each point (/, b) on the sky, one can determine 2 functions, p(l, b) which is the photon- 
axion conversion probability integrated along that line of sight within the galaxy and A(l, b) 
which is the effective HiRes stereo exposure in the same direction. The latter direction- 
dependent exposure is determined as described in Ref. |23j, based on the distribution of 
HiRes Monte-Carlo events in the zenith and azimuth angles and sidereal time. 

We use the list of correlated BL Lacs from [23] and compare their distribution with respect 
to axion conversion probability with the full sample of 156 BL Lacs studied in [231 [21] ■ The 
correlated events are plotted, together with the exposure, in the left panel of Fig. El One 
sees that they do not appear to follow exactly the expected random distribution (they would 
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be more likely to turn up in the most densely shaded region if the distribution of sources was 
isotropic). On the other hand, in the right panel of Fig. [U]we plot the expected distribution 
for the axion scenario which is proportional to the product p*A and notice that the agreement 
at least appears to be better. 

To quantify these suspicions, we use the method recently becoming popular in tests of 
global anisotropy of UHECR arrival directions [101 HD H2] • For each BL Lac with coordinates 
(h,bi) we calculate f\ = p(li,bi)A(li,bi) and Ai = A(li,bi). Then we compare, by means of 
the Kolmogorov-Smirnov test, the distributions of these Ai for BL Lacs which are correlated 
with cosmic rays and for all BL Lacs, correlated or not. The test gives a probability of 0.024 
that the two distributions of Ai are realizations of the same distribution, thus disfavoring 
the idea that correlated events come from random/isotropic regions on the sky. On the 
other hand, if we compare the distributions of fi for the two data sets, the probability is 
0.168, which means the distributions are compatible. We use the full BL Lac sample instead 
of simulated isotropic samples to account for non-uniform distribution of BL Lacs in the 
catalog. 

We have shown therefore that there does indeed appear to be a correlation between 
the arrival directions of neutral ultra high energy cosmic rays and regions of high axion- 
photon mixing probability. Before making any conclusions it is possible to list some other 
implications if such an effect turns out to be true. 

7 Consequences and future tests 

We have shown in the previous section that there is some motivation for a possible inter- 
pretation of the neutral events correlated with the position of BL Lac objects in the sky 
being due photons that have been able to traverse the Universe due to their conversion into 
axions in the source and back into photons in the Milky Way. More data with regard to the 
magnetic field in the Milky Way and more cosmic ray events will be able to add or subtract 
confidence in this interpretation but the scenario leads to several other consequences which 
may be tested in future studies. 

Primary particle type of the correlated events. Clearly, if axions are the explanation, 
then the primary particles of the correlated events should be photons. Currently studies 
of the primary particle type for the HiRes events are not published. The photon-primary 
hypothesis agrees perfectly with the absence of correlations in the Auger surface-detector 
data [13] : the photon energies are underestimated [28] by this detector by a factor of four on 
average [29J , so that the correlated events would be lost among a large number of hadronic 
events of lower energy. Such a situation should also be the case in the future data, although 
as alluded to earlier, the correlations should be seen in the data of fluorescent detectors of 
Pierre Auger and Telescope Array and in the surface detector of Telescope Array. 
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Secondary photons and the extended image. While 1/3 of the original gamma-ray flux 
is converted to axions in the source, the remaining 2/3 of photons rapidly lose their energy 
in electromagnetic cascades on the background radiation. As can be seen from Fig. HI this 
happens on distances d of order several Megaparsecs, that is, for distances D to BL Lacs 
of order hundreds Megaparsec an extend gamma-ray halo should form for EGRET/Fermi 
energies. The angular size of the halo is then 6 « which is a fraction of a degree, 
well below the width of the point spread function of either EGRET or Fermi. This means 
that the correlated BL Lacs should be gamma-ray sources in agreement with the results 
of Refs. |m H5]. The method of detection of the size of extended images of this kind is 
discussed in Ref. |4"B"] . 

The flux of a single source may be determined as follows:- The number of sources on 
the sky determined by the statistics of clustering [UJ is N s ~ 60. The correlated events 
constitute roughly 0.03 of the total number of events with E > 10 19 eV [48J. Given the 
total HiRes-detected flux at 10 19 eV g9] of JE 3 « 2 x 10 24 eV^m-^sr -1 , we estimate the 
single-source flux at 100 MeV as ~ 4 x 10~ 9 cm _2 s _1 , well within the Fermi sensitivity. 

Axion parameters. The model requires the axion-like particle with mass m < 10~ 7 eV 
and the inverse coupling to photon close to the current experimental limits, M ~ (1 -j- 10) x 
10 10 GeV. The most direct confirmation of the scenario would come from the discovery of 
that particle. This region of the parameter space is available for exploration with CAST at 
sufficiently large exposure. The axion with these parameters may also affect the polarization 
of extragalactic radio sources [501 EI] ■ 

8 Conclusions 

The existence of an Axion like particle with an inverse coupling M ~ 10 10 GeV and a low 
mass m < 10~ 7 eV has been invoked by other authors to explain the detection on earth of 
TeV photons from cosmological sources - flux which is difficult to explain given that such 
photons should produce electron-positron pairs on the cosmic infrared background [12] . The 
presence of such a particle would enable some photons to convert into axions in the source, 
travel over the intervening space without interacting with the background radiation before 
turning back into photons in the magnetosphere of the galaxy. Those authors showed that 
for a particular field profile, namely the HMR profile described in section 0of this paper, the 
presence of such a particle coupling to photons could explain how such high energy photons 
could reach earth. 

In this work we have tried to use the same method to explain a set of ultra high energy 
cosmic ray events which seem to come from BL Lac objects. Since such events seem to 
lead straight back to the source, the particles should be neutral because charged particles 
would be deflected by the magnetic field of the galaxy. However, no neutral particles in the 
standard model seem capable of traversing the Universe at such a high energy. The same 
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idea of photons turning into axions in the source and back into photons in the galaxy is 
immediately applicable to this cosmic ray situation. 

We test this hypothesis by assuming the same axion parameters and magnetic field as 
those that were used to successfully explain the detection of TeV photons from blazars 
in [12] . We use a statistical test which takes into account the varied sky coverage of the 
cosmic ray experiments to find out if the location of the sources of apparently neutral cosmic 
rays coincides with the directions in the sky where the probability of axion-photon mixing is 
highest. We do find a positive correlation and estimate that the probability of this correlation 
occurring through chance is only 2.4%. 

Clearly, more data is needed to show whether or not this correlation has occurred by 
chance, both cosmic ray data and data on the arrival of TeV photons from cosmological 
sources will help to add support to or rule out this hypothesis. 

It has been suggested that the existence of such an axion like particle would also help ex- 
plain some other astrophysical conundrums such as the white dwarf luminosity function [52] 
and large-scale correlations which may exist in the polarization of quasars, e.g. |3T]. 

Finally it is appropriate to re-iterate that the parameters of interest for this effect suggest 
a weak coupling for the axion, but not so weak that it cannot be probed by experiments such 
as CAST. The low mass required to ensure that one is in the region of maximal mixing means 
that such an axion should be able to be ruled out or confirmed by CAST as it continues to 
take data. 

We are grateful for discussions with V. Rubakov and I. Tkachev. This work was supported 
in part by the grants RFBR 07-02-00820 and NS-1616.2008.2 (ST). 
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